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Abstract. In the present article, we simulate the propagation of binary sig- 
nals in semi-infinite, mechanical chains of coupled oscillators harmonically 
driven at the end, by making use of the recently discovered process of nonlinear 
supratransmission. Our numerical results — which are based on a brand-new 
computational technique with energy-invariant properties — show an efficient 
and reliable transmission of information. 



The process of nonlinear supratransmission consists of a sudden increase in the 
amplitude of wave signals transmitted into a nonlinear chain by a harmonic dis- 
turbance at the end, irradiating at a frequency in the forbidden band gap. The 
phenomenon was first discovered in mechanical chains of oscillators described by 
coupled sine- Gordon and Klein- Gordon equations pQ, and it was quickly studied in 
other nonlinear models [2J |3l |4j [5] . Several applications of nonlinear supratransmis- 
sion have been realized so far [6l [7] , and more physical applications have been given 
a mathematical foundation to its realization [8l [9] have been suggested; however, 
the problem of transmitting coded information into nonlinear chain systems has 
not been attacked at all. 

In this paper, we develop an application of nonlinear supratransmission to the 
propagation of binary signals in weakly damped semi-infinite mechanical chains of 
coupled oscillators, by modulating the amplitude of the driving signal at the end. 
The first section introduces the mathematical model under study and the energy 
expressions to be used. A numerical study of the idealized problem is carried 
out in the following section; here we provide a numerical analysis of the breather 
propagation, based on the method proposed in [10j . which in turn is based on a 
numerical scheme to approximate radially symmetric solutions of modified Klein- 
Gordon equations . The next section presents a simulation of the propagation 
of a binary signal in a weakly damped system. Finally, we present a section of 
concluding remarks and propose complementary directions of research. 

1. Preliminaries 

1.1. Mathematical model. Throughout this paper, we assume that a, f3 and 7 
are nonnegative real numbers, and that c > 1 (see pQ). Likewise, we consider a 
system {unj^i of oscillators satisfying the mixed- value problem studied in |10) . 
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Figure 1. Bifurcation diagram of critical amplitude vs. driving 
frequency for system ([T]) for various values of 7: (solid), 0.1 
(dashed), 0.2 (dash-dotted), 0.3 (dotted). 



namely, 



(1) 



dt 2 
subject to : 



"dt 



2 . a du n 
x P~dt 
«n(0) = 0, 



V'(u n ) = 
n e Z+, 
n e Z+, 
t > 0, 



0, 



where c is the coupling coefficient, and a and j3 evidently play the roles of internal 
and external damping coefficients, respectively. Here, A^u„ is used to denote the 
spatial second-difference it n +i — 2u n + u n -i for every n £ Z + , the boundary-driving 
function is given by ip(t) = A(t) sin(f2t) for every t G (0,+oc), and V(u n ) 
1 — cos(m„) — 7it„ where, due to the analogy with the Josephson model [12| . 7 will 
be called here the normalized current. Notice that the Hamiltonian of the rt-th 
lattice site is given by 



Hm — — 

2 



c 2 (u, i+ i - u n ) 2 } + V(u n ), 



for any differentiable function V. After including the potential energy from the 
coupling between the first two oscillators, the total energy of the system becomes 



1.2. Numerical schemes. We consider a finite system of N differential equations 
satisfying |T]), and a regular partition — to < t\ < ■ ■ ■ < fjvf = T of the time 
interval [0, T] with time step equal to At. For each k — 0,1, ... , M, let us represent 
the approximate solution to our problem on the n-th lattice site at time i/. by u^- 



If we convey that <5 t it„ = 

^< = <+i ~ 2< + in- 
discrete form 



,fe+i 



,fc-i 



that Sivt 



,fc+i 



2«* 



,fc-i 



and that 



j , the differential equations in our problem take then the 
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Figure 2. Local energies H n (t) of (JT]) vs. lattice site n and time 
t, corresponding to the transmission of binary signal '1111'. 



(2) 
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2At 



.fc+1 



k-1 



= o, 



where (3' includes both the effect of external damping and a simulation of an ab- 
sorbing boundary slowly increasing in magnitude on the last N — No oscillators. 
More concretely, we let UN+i(t) be equal to zero at all time t, and let f3' be the 
sum of external damping and the function 



P"(n) = 0.5 



1 + tanh 



2n - N + N 
6 



where usually iV = 50 and N > 200. 

Finite-difference scheme ^fy (which is a modified version of the one proposed in 
[11] ) is consistent with our mixed- value problem and conditionally stable, having 
the inequality (cAt) < 1 + (a + /?' /4) Arj as a necessary condition for stability 
when /3' is assumed constant [10 . Moreover, if the energy of the system at the k-ih 
time step is computed using the expression 



1 M fn k 



At 



2 M 



A I 



then the discrete rate of change of energy turns out to be a consistent approximation 
of order O(At) 2 for the corresponding instantaneous rate of change. 



2. Numerical study 



4 



J. E. MACIAS-DIAZ AND A. PURI 



500 




J 00 - 



450 - 



/ 



350 - 



£ 300 - 



? 250 - 



m 200 - 



150 - 



/ 



1 DO - 



50 - 



100 



200 



300 
Time 



400 



500 



600 



Figure 3. Time-dependent graphs of the position of the breathers 
generated in (pj by the binary signal '1111'. 

2.1. Bifurcation analysis. The existence of a bifurcation threshold of the en- 
ergy administered into a semi-infinite chain of damped coupled oscillators de- 
scribed by ([T]) has been established and numerically predicted in [10] for a po- 
tential V(u) = 1 — cosu. Numerical experiments on undamped mechanical chains 
with nonzero normalized bias currents have shown that the process of nonlinear 
supratransmission is likewise present in these models. In fact, Fig. [T] provides 
bifurcation diagrams of driving amplitude at which supratransmission first occurs 
vs. driving frequency Q for a system of 200 undamped oscillators with coupling 
coefficient equal to 4 and several values of 7, over a time interval [0,T(f2)] where 
T(J7) is equal to 200 for all frequencies except for those satisfying il > 0.95, in 
which case T"(fi) had to be increased up to 500. 

2.2. Moving breather solutions. For the sake of simplification, let 7 = and 
consider a discrete system described by ([T]), harmonically driven at the boundary 
by a frequency SI. A binary bit b will be transmitted into the medium during a 
fixed and sufficiently long period of signal generation P equal to an integer multiple 
of the driving period, by defining 



for every t £ [0, P], where A s represents the critical amplitude at which nonlinear 
supratransmission starts and C > is an adjusting constant. 

Fig. [U shows the time evolution of the local energies H n corresponding to so- 
lutions of the undamped problem JT]) with coupling coefficient equal to 4, driving 
frequency 0.9 (in which case A s is approximately equal to 1.79) and C = A s . The 
period of signal generation is equal to 20 driving periods, and the binary code trans- 
mitted at the boundary is '1111'. The results show that a single moving breather 
is generated per period, and that the phase velocity v p through the medium is 
approximately constant. To verify this claim, we include in Fig. [3] the graphs of 
the positions of the four breathers generated by the binary signal vs. time. The 
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Figure 4. Time-dependent graphs of the maximum energy of the 
breathers generated in ([1} by the binary signal '1111'. 

fact that the breathers are transmitted at a constant velocity is now obvious, the 
common phase velocity being approximately 2.549. 

In order to determine the strength of the emitted signals at sites located far away 
from the source, it is important to determine the time behavior of the maximum 
local energy attained by a discrete breather. Fig. |4] presents this behavior for the 
breather solutions obtained under the conditions above. Notice that the maximum 
local energy is well above a cutoff limit of 2. 

3. Application 

The discrete medium described by ([1} with a~(3 = ^ = lx 10~ 3 , and the 
rest of the parameters as in the previous section, will be our object of study in 
this section. Local energy-based reception devices will be placed on the 100-th and 
300-th lattice sites, and the signal to be transmitted through the chain system is 
'10111001011011101001'. Our system will consist of 600 sites, and a time step of 
0.05 will be employed. 

Assuming the general convention that site no > 100 will start signal reception 
at time t — (no — 100) /v p , Fig. [5] shows the evolution of the local energy in the 
100-th and 300-th sites in terms of the time normalized with respect to the period 
of signal generation. In either case, it is clear that the transmission of a bit equal to 
1 in the n-th period is completely characterized in the graph by peak(s) of height 
greater than the cutoff limit 2 on the interval [n — 1, n] . 

4. Conclusions and perspectives 

In this letter, we have proposed a simple mathematical model to transmit bi- 
nary information in discrete, semi-infinite chains of couped oscillators using the 
process of nonlinear supratransmission. In the absence of dispersive and dissipative 
effects, our model (which is based on the modulation of amplitudes of source sig- 
nals with constant frequency) has shown to be highly reliable for sufficiently long 
periods of single-bit generation, independently of the distance between the source 
of transmission and the point of reception. 
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Figure 5. Local energy of the 100-th (left) and 300-th (right) sites 
in (fTJ) vs. normalized time, as a response to the transmission of 
the binary signal '10111001011011101001'. 

When weak damping is present the general picture does not change much. 
Stronger damping, however, manifests itself through a substantial decrease in the 
amplitude of the maximum local energy of the moving breathers with respect to 
the lattice site. In a forthcoming work, we will examine the possibility to overcome 
this problem via the concatenation of chain systems, where the driving at the be- 
ginning of each lattice will irradiate with an amplitude equal to a value just below 
its critical point multiplied by the amplitude of the last site in the previous lattice. 

Finally, we wish to point out that the problem of determining whether it is 
possible to design a propagation system of binary signals in Josephson junction 
arrays still remains an open topic of research. Moreover, in view of the recently 
discovered phenomenon of nonlinear infratransmission (or lower-transmission, as 
named by the authors [3]), the problem of finding more efficient pathways to achieve 
signal transmission in other models is still an open question of general interest. 
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